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This paper presents the results for the charged-particle multiplicity and transverse 

energy distributions at mid-rapidity in Au - Au collisions at y/s NN = 130 GeV measured 

with the PHENIX detector at RHIC. The values of dN c h/dr]^ =0 and dET/dr]\ v=0 , analyzed 

as a function of centrality, show a consistent steady rise. For the 5% most central collisions 

f^ they are ~ 70% larger compared to the SPS results for Pb-Pb collisions at y/s NN = 

17.2 GeV. The ratio of (Er r )/(N c f l ) remains constant as a function of centrality at 0.8 GeV, 

as also observed at the SPS at CERN and at the AGS at BNL. 

^> 

V ■ 1. INTRODUCTION 

Global variables such as dN c h/dr]\, n=0 and dET/drj\ v=0 are important for the charac- 
terization of high-energy nuclear collisions. They constrain theoretical models, help to 
discriminate among various mechanisms of entropy and particle production and allow 
to derive information about the initial conditions [|TJ. During the first run of RHIC at 
y/s NN = 130 GeV charged particle multiplicity and transverse energy density distributions 
were measured using the PHENIX detector. The results are analyzed as a function of 
centrality, defined by the number of nucleons participating in the collision and compared 
to similar results obtained at the much lower energies of the CERN SPS and the BNL 
AGS. 

2. THE PHENIX DETECTOR 

A detailed description of the PHENIX detector can be found in f2|. The detector 
consists of an axial-field magnet with two central arms, the East and the West arms, 
each one subtending 90° in azimuth (0) and ±0.35 units of pseudo-rapidity {rf). Two 
spectrometers, the North and South spectrometers, at forward angles (10° < 9 < 35°) 
around the beam axis serve to identify and track muons. 

The present analysis relies primarily on four PHENIX subsystems: the Zero Degree 
Calorimeters (ZDC) and the Beam-Beam Counters (BBC), used to derive the trigger and 
the off-line event selection, two layers of Pad Chambers (PC), called PCI and PC3 in 
the East arm, used to determine the charged particle multiplicity and two sectors of the 
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Electromagnetic Calorimeter (EMCal) in the West arm used to measure the transverse 
energy. These detector subsystems are shown in Fig.|l|. 



PCI 

East West 



EMCal 




Figure 1: PHENIX detector subsystems used in this work. 
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The ZDC are two small 
transverse area hadron calorime- 
ters that measure the neutron 
energy within \rj\ > 6. Be- 
sides their role as triggering 
devices they contribute to the 
centrality measurement, be- 
cause their energy signal is 
proportional to the number of 
spectator neutrons ||. The 
BBC comprise of two arrays of 
64 quartz-radiator Cherenkov 
detectors, located at ±1.44m 
from the center of the interac- 
< 3.9. The BBC also provide the start timing 
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tion region, covering 2tt in 
and vertex position ||. 

The PC provide three-dimensional coordinates along the charged-particle trajecto- 
ries 0. Each PC layer has 8 wire chambers with cathode pad readout. PCI and PC3 
are placed at radial distances of 2.49 m and 4.98 m, respectively. The instrumented part 
of the EMCal is a Lead/Scintillating tile Calorimeter || which covers the pseudorapidity 
range \rj\ < 0.38 with an azimuthal aperture of A0 = 45°. The depth of the calorimeter 
corresponds to 18 radiation lengths, equivalent to 0.85 interaction lengths for hadrons. 

The event selection consists of an interaction trigger generated by a coincidence be- 
tween the north and south BBC with at least two detectors firing in each of them and a 
requirement of the collision vertex position, to be within \z\ < 17 cm. Based on detailed 
simulations of the BBC, this trigger reflects [92 ± 2(syst)]% of the nuclear interaction 
cross section of 7.2 barns. More information on the trigger selection can be found in [[?]]. 
The data sample used in this analysis consists of 137,784 events taken with zero magnetic 
field. 



3. ANALYSIS 

The number of primary charged particles per event was determined on a statistical 
basis by correlating hits in PCI and PC3 in the East arm, rather than by explicit track 
reconstruction, using the following algorithm: all hits in PC3 are combined with all 
hits in PCI and the resulting lines are projected onto a plane through the beam line, 
symmetrically placed with respect to the tracking arm. This distribution of the distance 
R from the track intersection point to the event vertex is shown in Fig.^. Two main 
classes of tracks contribute to it: real tracks (peaking at low R) and tracks from the 
obvious combinatorial background inherent to the adopted procedure. The latter can be 
determined by a mixed event technique; in the present analysis, each sector in PCI was 
exchanged with its neighbor and the resulting combinatorial background is shown in Fig.||| 
by the dotted line. The yield of this background increases quadratically with R (leading 
to a linear dependence in the differential dN/dR vs. R presentation of Fig.^j). 



The R distribution of real tracks is obtained by subtracting an average background on an 
event-by-event basis from the total number of tracks. The sharp peak in the distribution 

at small R is due to tracks from primary 




background 



real 



10 



15 



20 



25 
R [cm] 

Figure 2: Track distance from vertex. 



iacK distance uom veitex particles originating at the vertex and the 

long tail to in-flight decay products of pri- 
mary particles. 

The total number of real tracks deter- 
mined after background subtraction has to 
be corrected for several factors: corrections 
of 15.3% and 1.2% account for inactive re- 
gions and inefficiency in both PC layers, 
respectively. The number of tracks outside 
the acceptance window (25 cm) amount to 
4.3%. The track losses due to the finite 
double hit resolution of the PC depend on the event multiplicity. The losses are due 
to two effects: the direct counting of tracks (13.3% for the 5% most central collisions) 
and the subtraction of the combinatorial background by event mixing (3.6% of the back- 
ground). Finally, a 2.8% correction takes into account the effect of charged and neutral 
particle decay in flight. A geometrical scaling factor of 5.82 allows to convert the corrected 
number of tracks to the charged particle multiplicity in one unit of pseudorapidity. The 
systematic errors were thoroughly studied using both the data and the guidance from 
detailed Monte Carlo simulations. The total systematic error amounts to 6.5% at the 
highest multiplicities. See for more details about the analysis. 

The absolute scale of the EMCal was pre-calibrated during construction by cosmic- 
ray muons and in-situ by monitoring the minimum ionizing peak of relativistic charged 
hadrons and identified high momentum electrons. The ir° mass, reconstructed from pairs 
of EMCal clusters of total energy greater than 2 GeV is shown in Fig.[| and is found to 
be within 1.5% of the true value. 

The transverse energy E? was computed 
for each event using clusters of minimum 
energy of 20MeV, from individual towers 
with deposited energy of more than 3MeV. 
The angle of the cluster was computed 
from the centroid of its energy assuming a 
particle originating from the event vertex. 
A correction factor of 1.17 converts the 
measured E T to the total E T in the fidu- 
cial aperture of the EMCal, and includes 
the response of the detector to charged and 
neutral particles from the event vertex and 
for energy in- and out-flow from the fidu- 
cial aperture. This factor was determined from a full Monte Carlo simulation of the 
detector, using HIJING || as event generator. A correction of 3% accounts for disabled 
calorimeter towers. A scaling factor of 10.6 converts the total measured Et value into the 
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Figure 3: ir° peak reconstructed by EMCal. 



transverse energy in one unit of pheudorapidity. The systematic error of the measurement 
was determined to be 4.5% independent of multiplicity. For more details see 0. 

Centrality classes were determined based on the ZDC vs. BBC response and related 
to the number of participating nucleons (N p ) and binary collisions (N c ) using a Glauber 
model 0. 

4. RESULTS 

The measured distributions of dN c h/drj and dEx/dr] at mid-rapidity are shown in Fig.f|. 
Both quantities show very similar behavior. The peak at low values corresponds to pe- 
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Figure 4. Multiplicity (left) and Transverse energy (right) densities at mid-rapidity 



ripheral events, the middle region to mid-central collisions and the fall-off to the most 
central events. The shape of the fall-off is mainly determined by the limited detector 
acceptance rather than by real fluctuations. The lower axis in each plot corresponds to 
the measured distribution calculated for one unit of pseudorapidity. The four bell-shape 
curves in each plot show the distribution of the 20% most central events in steps of 5%. 
Fig.|5] shows the multiplicity and the transverse energy per pair of N p vs N p . One can 
see that the extrapolation to N p = 2 of the measured dN c f l /dn ri= Q points is in agreement 



with the UA5 data for pp flTOf. The figure also shows the results obtained by WA98 ITT 
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Figure 5. dNch/dr)\ v= o (left) and dE? / 'dr)\ v =o 

a ^ \ /s nn = 17.2 GeV. The curves represent a fit to the data points with dX/dr] oc N p 
(X = N ch or Et) as used in fTTj. For dN ch /drj one finds a = 1.16 ± 0.04 for PHENIX 



whereas a = 1.07±0.04 and a = 1.05 ±0.05 for WA98 and WA97/NA57 pf respectively. 
For dE T /dT] the PHENIX value is a = 1.13 ±0.05 compared to a = 1.08 ±0.06 for WA98. 
One can conclude that at RHIC energies both measured quantities show a consistent 
increase with centrality, with an a value clearly larger then one, whereas at the SPS 
energy a was found to be very close to 1. 

Fig.|5] compares the measured dN c h/drj 

with the theoretical predictions of the 

EKRT [13 1 saturation model and HI- 
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Figure 6: Comparison with models. 
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N p is in contrast to the prediction of the 
EKRT model, and in qualitative agreement 
to HIJING, however the latter underpre- 
dicts the yield by ~15%. Following HI- 
JING, one can assume that the increase is 
due to "hard" processes that scale with N c 
and fit the results with the parameteriza- 
tion A ■ N p + B ■ N c . One finds that the 
ratio B/A for dN ch /dr] is 0.38 ± 0.19. For 
the dEx/dr) the ratio B/A = 0.29 ± 0.18. The large errors are due to the negative corre- 
lation between A and B. Therefore, in the framework of the HIJING model the "hard" 
contribution grows from ~30% for the 45-50% centrality bin to ~50% for the most central 
bin (0-5%). 

Rapidity densities per N p obtained in different experiments are shown in Fig.|7| as a 
function of the center of mass energy. The compiled data were measured for different 
centrality percentiles (values shown in brackets) in the center of mass system for PHENIX 
and PHOBOS or laboratory system for all other experiments. In making this figure it was 
assumed that dX/dy ~ dX/dr) in the lab system and factors of 1.19 at y/s NN = 130 GeV 
and 1.24 at y/s NN = 56 GeV were used to account for the r] — > y transformation in 
the cm. system. The PHENIX result is in good agreement with the recent PHOBOS 
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7. dN ch /dy\ y= o (left) and dE T /dy\ y=0 (right) per pair of N p vs.^/s 
Data are taken from: PHENIX 0,|, PHOBOS |E5 
0,13, E814/E887 PH, and E802 [g(]]. 
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measurement. Comparing the mean values for RHIC and SPS shows that the particle 
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and transverse energy production per participant increase by ~ 70%. The data points 
are consistent with dX/dy oc ln(-Js NN ) behaviour over a broad range of energies. 

Finally, Fig.[|] shows the ratio (Et) / (N c h) of the transverse energy per charged particle. 
The ratio stays constant, equal to 0.8 GeV over a broad range of centralities, as was also 
observed by WA98. Moreover, the ratio, stays also constant as a function of y/s NN , from 
AGS to SPS and up to RHIC energies. 
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(Et) I (Nch) vs. N p (left) and vs. y/s NN (right) for 0-5% centrality bin. Data 
from: PHENIX §g, WA98 Q, NA49 0,0], and E814/E887 MM- 
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